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ABSTRACT 

Using  selectively  detected  laser-induced  fluorescence,  the 
rotational  state  distribution  of  the  Bal  product  has  been  measured  for 
the  beam-gas  reaction,  Ba  +  HI  -*  Bal{D=8)  +  H.  Owing  to  the  highly 
constrained  kinematics  for  this  system,  these  measurements  can  be  used 
to  derive  the  reaction  probability  as  a  function  of  the  impact  parameter 
for  this  channel,  called  the  “specific"  opacity  function,  once  the  reaction 
probability  as  a  function  of  velocity  has  been  determined.  Unfortunately, 
lack  of  knowledge  of  the  exoergicity.  and  the  height  of  any  energy 
barrier  prevents  a  conclusive  determination  of  the  specific  opacity 
function  for  this  reaction.  Instead,  various  approximate  opacity 
functions  are  estimated  based  on  different  models  of  the  velocity 
dependence  of  the  reaction  channel  studied.  If  the  reaction  probability  is 
the  same  for  all  relative  collision  velocities,  then  the  Bal(s=8)  specific 
opacity  function  peaks  strongly  near  2.6  A  with  a  full  width  at  half 
maximum  of  1.0  A.  However,  the  possible  presence  cf  a  small  energy 
barrier  in  the  entrance  channel  causes  a  cutoff  in  the  relative  collision 
velocity  distribution,  and  this  type  of  velocity  dependence  would 
significantly  affect  the  shape  of  the  specific  opacity  function. 
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_>tising  selectively  detected  laser-induced  fluorescence,  the  rotatlbrpl  state  distribution  of  the  Bal 
product  has  been  measured  for  the  beam-gas  reaction  Ba  +  HtiJ<'BaI(v  =  8)  -f-  H.  Owing  to 
the  highly  constrained  kinematics  for  this  system,  these  measurements  can  be  used  to  derive  the 
reaction  probability  as  a  function  of  the  impact  parameter  for  this  channel,  called  the‘~Specific*s- 
opacity  function,  once  the  reaction  probability  as  a  function  of  velocity  has  been  determined. 
Unfortunately,  lack  of  knowledge  of  the  exoergicity  and  the  height  of  any  energy  barrier  prevents 
a  conclusive  determination  of  the  specific  opacity  function  for  this  reaction.  Instead,  various 
approximate  opacity  functions  are  estimated  based  on  different  models  of  the  velocity  dependence 
of  the  reaction  channel  studied.  If  the  reaction  probability  is  the  same  for  all  relative  collision 
velocities,  then  the  Bal(  v  =  8)  specific  opacity  function  peaks  strongly  near  2.6  Awith  a  full 
width  at  half-maximum  of  l.CUVHowever,  the  possible  presence  of  a  small  energy  barrier  in  the 
entrance  channel  causes  a  cutoff  in  the  relative  collision  velocity  distribution,  and  this  type  of 
velocity  dependence  would  significantly  affect  the  shape  of  the  specific  opacity  function. 

A 


I.  INTRODUCTION 

Consider  the  elementary  reaction 

A  4-  BC  — *  AB(v)  +  C,  (1) 

where  A,  B,  and  C  are  atoms  and  v  denotes  the  vibrational 
state.  Although  the  AB  product  is  often  detected  indepen¬ 
dent  of  the  vibrational  state,  the  use  of  laser  spectroscopic 
techniques  allows  experimentalists  to  study  the  reaction 
channel  for  an  individual  product  vibrational  state.  Regard¬ 
less  of  whether  a  specific  vibrational  state  is  detected,  con¬ 
temporary  experimental  studies  of  reaction  dynamics  tend 
to  ignore  the  detailed  dependence  of  the  reaction  probability 
on  the  impact  parameter  b,  which  is  defined  as  the  distance 
of  closest  approach  of  the  incoming  reagents  if  they  move  in 
undeflected  straight  lines.  This  is  primarily  a  consequence 
of  the  fact  that  experimentalists  cannot  control  the  choice  of 
the  impact  parameter  in  a  collision. 

The  impact  parameter  is  related  to  the  cross  section, 
a ,  ( ),  for  reaction  ( 1 )  by 

<7,.(t’„.| )  =  2tt  I  P,.(b,unl  )b  db,  (2) 

J() 

where  Pv  (b,vrel )  gives  the  reaction  probability  for  a  given 
and  b  to  produce  the  AB  product  in  a  specific  vibrational 
level,  v  (averaged  over  all  reagent  orientations).  Thus, 
P,(b,  vnl)  is  called  a  "specific”  opacity  function,  which 
must  be  distinguished  from  the  “overall”  opacity  function 
P(b,V"i ).  This  overall  opacity  function  describes  the  reac¬ 
tion  probability  for  a  given  i>„,  and  b  to  produce  the  AB 
product  in  any  vibrational  state,  and  is  related  to  the  specific 
opacity  functions  by 
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P(b,vrci)  =  ^PAb,  t»rc|).  (3) 

The  “velocity-averaged”  specific  opacity  function  is 
then  defined  as 

PAb)  =  f  P,  ( b,vrcl )/( i’rd  )dvrvi /  f  /( i>r,,|  )dvKl ,  (4) 

Jo  Jo 

where /(y„| )  is  the  relative  velocity  distribution  of  the  rea¬ 
gents  in  the  center  of  mass  collision  frame.  In  terms  of  P,  ( b ) , 
the  velocity-averaged  cross  section  a,,  is  related  to  the  specif¬ 
ic  opacity  function  by 

a,  =2tt\  P,\b)bdb.  (5) 

J() 

As  first  pointed  out  by  Herschbach,'  under  certain  fa¬ 
vorable  circumstances  it  is  possible  to  remove  the  seemingly 
intrinsic  average  over  b  and  glimpse  the  detailed  role  of  the 
impact  parameter  in  governing  the  outcome  of  a  reactive 
encounter.  The  key  is  to  choose  a  kinematically  constrained 
reaction,  for  example,  that  of  a  heavy  atom  plus  a  heavy- 
light  diatomic  reacting  to  give  a  heavy-heavy  diatomic  pro¬ 
duct  plus  a  light  escaping  atom. M  A  near  limiting  example  of 
such  a  reaction  is 

Ba -f  HI  — *  Bal(v)  +  H.  (6) 

As  a  consequence  of  the  conservation  of  angular  momen¬ 
tum,  nearly  all  of  the  reagent  angular  momentum  L  is  chan¬ 
neled  into  rotational  angular  momentum  J  of  the  diatomic 
product.  Here  |L|  =  pvrclb,  where  for  reaction  (6), 

\/p  =  l/mB>  +  1  /m,u.  (7) 

Thus, /or  a  given  vn.,,  the  distribution  of  the  reactive  impact 
parameters  is  mapped  into  the  population  of  the  product 
rotational  levels  through  the  relation 

I  J)  =  |L|  =pvrc,b.  (8) 

Although  this  possibility  of  determining  the  specific 
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I 


opacity  function  for  kinematically  constrained  systems  was 
proposed  more  than  20  years  ago,  ’  there  has  been  no  attempt 
to  carry  out  the  actual  measurements  to  date.  This  is  mainly 
due  to  the  fact  that  the  reaction  product  must  be,  by  defini¬ 
tion,  a  heavy  molecule,  and  consequently,  the  spectra  are 
very  congested,  hindering  a  rotational  analysis.  We  present 
here  the  first  experimental  measurement  of  the  product  rota¬ 
tional  distribution  for  the  Ba  +  HI  — ►  Bal(v  =  8)  +  H  sys¬ 
tem,  and  use  this  information  to  estimate  the  specific  opacity 
function  for  this  reaction  channel,  assuming  various  models 
for  the  dependence  of  the  reaction  probability  on  relative 
collision  velocity. 

When  the  relative  velocity  distribution  has  some  spread, 
the  product  rotational  population  is  affected  by  the  velocity 
dependence  of  the  reaction  probability.  This  velocity  depen¬ 
dence  can  be  separated  into  (a)  the  probability  that  a  trajec¬ 
tory  reaches  a  transition  state,  and  (b)  the  probability  that  a 
reaction  takes  place  after  a  trajectory  has  reached  a  transi¬ 
tion  state.  We  call  the  former  a  "transmission"  velocity  de¬ 
pendence,  and  the  latter  an  "intrinsic”  velocity  dependence. 
The  transmission  velocity  dependence  is  caused  by  the  exis¬ 
tence  of  an  activation  energy  barrier  in  the  entrance  channel 
and  by  the  energetics  of  the  reaction,  both  of  which  permit 
only  velocities  greater  than  some  threshold  to  react  to  pro¬ 
duce  the  Bal  product  in  v  =  8.  The  intrinsic  velocity  depen¬ 
dence  arises  from,  for  instance,  a  decrease  in  the  interaction 
time  between  the  reagents  at  high  collision  velocities  (see  pp. 
116-118  of  Ref.  2). 

As  will  be  discussed  in  Sec.  Ill,  the  lack  of  information 
on  the  exoergicity,  the  height  of  the  energy  barrier,  and  the 
intrinsic  velocity  dependence  of  the  reaction  probability  pro¬ 
hibits  a  unique  determination  of  the  specific  opacity  function 
for  this  reaction  system.  This  implies  that  the  relative  distri¬ 
bution  in  product  rotation  P(J)  can  only  be  related  to  the 
specific  opacity  function  Pv(b)  with  assumed  forms  of  the 
velocity  dependence. 


II.  EXPERIMENTAL 

The  experimental  apparatus  used  in  this  work  has  been 
described  elsewhere, ,_7  and  the  experimental  techniques  em¬ 
ployed  here  are  explained  in  detail  in  Ref.  7.  Thus  only  a  brief 
account  is  provided  here.  The  vacuum  chamber  consists  of 
two  differentially  pumped  chambers:  one  for  the  Ba  beam 
source,  the  other  for  the  beam-gas  reaction.  The  pressure  is 
measured  using  uncalibrated  Bayard-Alpert  type  ionization 
gauges.  A  typical  base  pressure  is  1 X  10~6  Torr.  A  stainless 
steel  crucible  containing  Ba  metal  (purity  99.5%,  Alfa  Pro¬ 
ducts)  has  an  orifice  of  0.4  mm  diameter  and  is  radiatively 
heated  to  985  °C.  Before  each  experiment,  HI  gas  (purity 
98.0%,  Matheson)  was  cooled  to  liquid  N2  temperature  and 
pumped  to  remove  H2  which  is  always  present  due  to  the 
equilibrium  2HI«=rH2  +  I2.  During  the  experiments,  the 
temperature  of  the  HI  cylinder  is  kept  at  —  20  °C  to  con¬ 
dense  I2.  At  this  temperature,  the  HI  pressure  is  on  the  order 
of  one  atmosphere,11  and  the  HI  gas  is  introduced  into  the 
reaction  chamber  via  a  double-needle  precision  valve.  The 
HI  pressure  in  the  reaction  chamber  is  maintained  at 
2.0X  10  4  Torr  (uncalibrated  ionization  gauge  reading). 

A  single-mode  ring  dye  laser  (Spectra-Physics  380D)  is 
pumped  by  the  5 14. 5  nm  line  of  an  Ar+  laser  (Coherent  CR- 
18).  Rhodamine  560  dye  is  used  with  a /?H  =  10  buffer  solu¬ 
tion  to  prolong  dye  lifetime.  The  typical  output  power  is  200 
mW  with  a  4  W  pump  laser  power.  The  laser  beam  intersects 
the  Ba  beam  at  right  angles.  The  resulting  fluorescence  is 
monitored  by  a  cooled  photomultiplier  (Centronics 
Q4283RA,  extended  S-20  photocathode)  through  a  1  m 
monochromator  (Interactive  Technology)  equipped  with  a 
1 200  grooves/mm  grating,  and  measured  by  a  standard  pho¬ 
ton  counter-rate  meter  combination.  A  Y AG-pumped  dye 
laser  (Quanta-Ray  DCR-1A/PDL-1 )  was  used  to  obtain 
broadband  excitation  spectra,  as  shown  in  Fig.  1 .  The  spec¬ 
trum  shown  in  Fig.  1  is  similar  to  the  one  reported  by  Cruse, 


FIG.  I.  A  low-resolution  excitation  spec¬ 
trum  of  the  Bal  C3II  system, 

using  a  Y AG-pumped  dye  laser  with  a 
bandwidth  of  ~0.3  cm-'.  Since  the 
Franck-Condon  factors  for  Av  =  0  tran¬ 
sitions  are  almost  unity  and  all  the  transi¬ 
tions  are  strongly  saturated,  the  intensities 
of  the  transitions  give  rough  estimates  of 
the  product  vibrational  populations. 
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FIG  2  A  3  cm  1  segment  of  the  Bal  C‘Ilj  *  (8*8)  excitation  spec¬ 

trum  ( upper  trace)  and  the  selectively  detected  laser-induced  fluorescence 
spectrum  (middle  trace)  An  F  fluorescence  spectrum  and  1  GHz  etalon 
markers,  which  were  used  to  calibrate  wavelengths,  are  shown  in  the  lower 
two  traces. 


Dagdigian,  and  Zare.1'  who  first  studied  this  reaction  under 
crossed  jet  conditions. 

While  the  large  mass  of  the  heavy-heavy  product  is  re¬ 
quired  for  a  rigorous  kinematic  constraint  on  the  reaction 
Ba  +  HI  — ►  Bal  +  H,  it  also  causes  the  visible  spectra  of  the 
Bal  product  to  be  congested  and  hence  difficult  to  ana¬ 
lyze.'’ For  example,  only  bandheads  and  rotational  con¬ 
tours  can  be  resolved  in  the  Bal  product  excitation  spectrum 
observed  in  the  beam-gas  reaction  Ba  +  HI  using  a  pulsed 
dye  laser  with  —0.3  cm  '  bandwidth  (see  Fig.  1 ).  When  a 
single-mode  dye  laser  with  a  bandwidth  narrower  than  the 
natural  linewidth  of  the  Bal  C  ~I  I  state  is  used,  the  excitation 
spectrum  (upper  trace  of  Fig.  2)  is  too  congested  to  permit 
an  accurate  measurement  of  the  rotational  distribution  or 


even  an  assignment  of  the  lines  in  many  cases.  The  small 
rotational  constants  of  the  upper  and  lower  states  of  the 
C2U-X21*  band  system  0.027  cm-1,  Ref.  6) 

combined  with  the  Rydberg  character  of  the  upper  state  that 
places  its  potential  almost  directly  above  and  nearly  parallel 
to  the  ground  state  potential  conspire  to  produce  this  con¬ 
gestion.  It  might  appear  then  that  the  unambiguous  mea¬ 
surement  of  a  product  state  distribution  for  this  reaction 
would  be  impossible.  We  have,  however,  succeeded  in  carry¬ 
ing  out  a  rotational  analysis6  7'1011  using  a  combination  of 
laser  techniques,  in  particular  optical-optical  double  reso¬ 
nance  and  selectively  detected  laser  induced  fluorescence 
(SDLIF).12  The  latter  method  is  of  particular  use  and  de¬ 
serves  discussion  here. 

Figure  3  shows  a  Fortrat  diagram  illustrating  the 
SDLIF  technique  for  the  Bal  C2n,/2-A’  21+  subband.  The 
six  possible  rotational  transitions  in  this  subband  can  be  sep¬ 
arated  into  two  groups  according  to  the  parity  component 
(A  doublet)  of  the  2II  excited  state.  Thus  excitation  via  the 
P,2  branch  results  in  fluorescence  in  the  Pl2,  (?,,  and  Rl2 
branches.  If  a  monochromator  isolates  the  fluorescence  of 
the  Q\  +  R,2  bandhead,  it  is  thus  possible  to  detect  selective¬ 
ly  excitation  of  individual  rotational  transitions  in  the  Pi2 
branch.  Figure  2  presents  a  3  cm  ' 1  segment  of  the  Bal  C-X 
excitation  spectrum  ( upper  trace )  and  the  selectively  detect¬ 
ed  fluorescence  spectrum  (middle  trace)  for  the  rotational 
levels  of  the  Bal(  v  =  8)  product  observed  in  the  Ba  -I-  HI 
beam-gas  reaction.  Note  that  the  rotational  distribution  is  so 
extensive  that  it  still  appears  flat,  even  when  over  35  rota¬ 
tional  levels  were  observed,  as  in  this  scan. 

Since  selective  detection  relies  on  isolation  of  fluores¬ 
cence,  it  is  mandatory  to  narrow  the  monochromator  slits, 
typically,  to  1 50//m.  This  unfortunately  causes  the  detection 
sensitivity  of  Pl2  excitation  to  vary  as  a  function  of  J,  since 
the  frequency  of  the  (?,  -t-  R]2  transition  changes  and  the 
"slit  function"  of  the  monochromator  is  not  constant  over 
the  Q ,  +  Ri:  bandhead.  This  problem  was  overcome  by  a 
direct  measurement  of  the  slit  function.  When  an  experiment 
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FIG.  3.  A  Fortral  diagram  of  the  Bal 
C *  system  Closed  dots  rep¬ 
resent  transitions  involving  /  level  upper 
states  (A -doubling  components >.  and  open 
dots  represent  transitions  involving  e  level 
uppc  states.  When  a  monochromator  col¬ 
lects  fluorescence  only  in  the  Q ,  +  Rtl 
bandhead  region,  the  excitation  of  the  P,; 
transitions  can  be  selectively  detected. 
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FIG.  4.  The  experimentally  observed  slit 
function  and  line  positions  of  the  Bal 
C-A'CS.S)  Qs  and  /?,.%  transitions. 


is  completed,  N,  gas  is  introduced  into  the  reaction  chamber 
and  scattered  light  from  the  laser  is  measured  as  the  laser 
scans  across  the  slit  function.  The  scattered  light  intensity 
divided  by  the  laser  power  gives  the  slit  function  of  the  parti¬ 
cular  experiment's  monochromator  setting  ( see  Fig.  4 ) .  The 
detection  efficiency  for  the  P]2U)  transition  DU)  is  then 
given  by 

D(J)  -  S  [QXU  —  1 )  ]  X  [slit  function  at  Q,(J  -  1 )  ] 

+  S[Rt,U  —  2 )  ]  X  [slit  function  at  R,;(7  —  2)  ], 

(9) 

where  5  denotes  the  line  strength  for  the  transition  in  paren¬ 
theses  divided  by  2J  +  1  and  the  slit  function  is  evaluated  at 
the  frequency  of  the  transition,  as  shown  in  Fig.  4.  Since  we 
are  interested  only  in  the  relative  intensities,  the  ratios  of  the 
line  strengths  are  used  in  Eq.  (9):  2  for  members  of  the  Q , 
branch,  and  1  for  members  of  the  /?,,  branch  (the  high-7 
limiting  ratios).1 1 

The  Bal  product  will  have  the  rotational  alignment,  de¬ 
fined  as  (P2(  J  •  v„i  )>14  where  J  •  vrd  is  the  cosine  of  the 
angle  between  the  product  rotational  angular  momentum  J 
and  the  reagent  relative  velocity  vector  vrel,  P2(x) 
=  ( 3  cos^  x  —  1 )  /2,  and  the  superscript  caret  denotes  a  unit 
vector.  Because  of  the  kinematic  constraint,  the  product  ro¬ 
tational  angular  momentum  is  equal  to  the  reagent  orbital 
angular  momentum  J  =  L  and  consequently  J  is  perpendic¬ 
ular  to  v„,  regardless  of  the  magnitude  of  the  angular  mo¬ 
mentum.  Hence,  the  rotational  alignment  is  independent  of/ 
and  the  effect  of  the  product  alignment  on  the  detection  sen¬ 
sitivity  of  each  rotational  state  can  be  negle  .  ted 

The  Bal  C  :II-.Y  21  *  band  system  is  known  to  be  satu¬ 
rated  very  easily.1' Ih  This  could  be  avoided  by  using  a  very 
low  laser  power  at  a  sacrifice  of  signal  intensity.  We  adopted, 
however,  a  different  approach.  The  unattenuated  laser  is 
used  to  nearly  saturate  the  transition  and  the  laser  power  is 
recorded  at  the  same  time  as  the  SDLIF  spectra  are  taken. 


Later,  the  dependence  of  the  signal  intensity  on  the  laser 
power  is  measured,  and  the  signals  are  corrected  accordingly 
(see  Fig.  5 1.  It  is  important  to  note  that,  since  the  laser  pow¬ 
er  is  very  high  and  the  transition  is  strongly  saturated,  the 
laser-induced  fluorescence  method  functions  in  this  instance 
neither  as  a  density  detector  nor  as  a  flux  detector,  but  as 
some  intermediate  between  the  two.  ^  Fortunately  the  corre¬ 
lation  between  t’rcl  and  ty  M  is  very  weak  in  the  Ba  -I-  HI 
beam-gas  reaction,  so  that  the  resulting  relative  distribu¬ 
tions  become,  within  the  experimental  error,  the  same 
whether  detected  as  a  flux  or  as  a  density. 

Recalling  that  |J|  =/r  it  is  important  to  determine 
i’rd  and,  therefore,  to  know  the  velocity  distributions  of  Ba 
and  HI.  The  HI  velocity  can  be  assumed  to  follow  the  Max- 
well-Boltzmann  distribution.  The  Ba  beam  velocity  is  mea¬ 
sured  directly  using  a  Doppler-shift  method.'7  which  em¬ 
ploys  two  laser  beams,  one  intersecting  the  Ba  beam  at  right 
angles  and  the  other  at  45°.  When  the  laser  is  scanned  across 
the  Ba  atomic  line  profile  ( [P°  -  'S„  at  18  060.264cm 
the  excitation  by  the  perpendicular  laser  gives  a  Doppler- 


Mti  5.  Laser  power  dependence  of  the  SDLIF  signal.  The  solid  s  tides  rep¬ 
resent  experimental  results. 
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Al*  (MH2) 


FIG  b  A  Dopplcr-shift  >ptctrum  ol'thc  Ba  'P  'S„  transition.  The  broken 
line  I'  the  ealeulateil  line  profile  using  the  Ba  velocity  distribution  gnen  by 
Eq  l  10). 


tree  signal,  while  the  excitation  by  the  laser  at  45°  prov  ides  a 
Doppler-shifted  signal,  as  shown  in  Fig.  6.  Note  that  the  Ba 
atomic  line  shows  several  components  due  to  different  iso¬ 
topes,  some  having  hyperfine  structure.1''  By  convoluting 
the  Doppler-free  line  shape  with  the  Doppler-shifted  line 
shape,  we  obtained  the  (unnormalized)  Ba  beam  velocity 
distribution  "  as 

/(t'n,  I  ■=  (I',,,  ):  e\p[  -  (t’„,  -  t\):/rr;],  ( 10 i 

where  t\  =  300  ms  1  and  a  =  (2kT/mtii> ) 1  12  at  T 
=  085  °C. 

III.  RESULTS  AND  DISCUSSION 

A.  Rotational  distribution  of  the  Bal(v  =  8)  product 

Figure  7  shows  the  rotational  distribution  of  the 
Ball  v  =  8)  product  produced  in  the  Ba  +  HI  reaction.  The 
highest  J  observed  here  is  350.5.  Thus,  in  order  to  probe 
more  than  300  rotational  levels,  the  single-mode  laser  must 
be  scanned  more  than  25  cm  '.  In  a  given  experiment,  each 
rotational  line  is  measured  at  least  twice,  and  the  average 
value  is  presented  here.  The  measurements  are  performed  at 
different  times  during  an  experiment,  in  an  attempt  to  avoid 
possible  systematic  deviations  due  to,  for  example,  a  long¬ 
term  variation  of  the  Ba  beam  flux.  The  error  bars  given  in 
Fig  7  are  associated  with  the  following:  ( 1 )  deviations  of  the 
measurements  for  each  J.  ( 2 )  deviations  of  measurements 
for  the  adjacent  ten  rotational  levels;  (3)  errors  from  the 
power  normalization:  and  (4)  errors  involved  in  the  slit 
function  corrections.  The  error  bars  should  be  taken  as  being 
on  the  safe  side. 

A  similar  experiment  was  performed  using  a  much 
weaker  laser  power  (  -  20  mW  )  to  test  the  effect  of  the  laser 
power  and  detection  scheme  ( flux  or  density )  on  the  result¬ 
ing  distribution.  The  results  of  this  experiment  agreed  with 
(lie  one  shown  in  F  ig.  7  within  experimental  error. 


B.  Specific  opacity  function:  No  velocity  dependence  ol 
reaction  probability 

From  Eq.  (8),  it  can  be  seen  that  the  rotational  distribu¬ 
tion  P(J)  provides  information  about  the  probability  of  re¬ 
action  as  a  function  of  both  vrcl  and  b.  The  distribution  of  vri!, 
=  Ivh»  —  v ui  ]  can  be  calculated  from  the  experimentally 
determined  f(vHj),  Eq.  (10),  and  a  Maxwell-Boltzmann 
distribution  for  the  HI  gas  using  a  Monte  Carlo  simulation 
technique.  The  simplest  treatment  is  to  ignore  any  depen¬ 
dence  of  the  reaction  probability  on  relative  collision  veloc¬ 
ity. 

In  order  to  determine  the  specific  opacity  function. 
P ,  (b),  from  the  experimentally  determined  P(J),  we  ex¬ 
pand  the  specific  opacity  function  Pv(b )  in  terms  of  8  func¬ 
tions,  as 

PAb)  =2.a,6(b  -  b,),  (II) 

where  b ,  is  some  impact  parameter,  and  a,  is  the  weighting 
factor  of  b,.  For  each  member  in  the  expansion,  the  corre¬ 
sponding  rotational  distribution  P,  ( J )  can  be  calculated  by 
using  Eq.  (8),  once  the  reactive  velocity  distribution  is  giv¬ 
en.  Thus,  the  observed  rotational  distribution  is  related  to 
the  calculated  distribution  by 

P(J)  =^a,P,(J)-  (12) 

Hence,  the  problem  of  obtaining  the  specific  opacity  func¬ 
tion  P,  (b)  is  reduced  to  determining  the  expansion  coeffi¬ 
cients  {<z,}  using  a  linear  least-squares  fitting  method. 

The  observed  and  generated  rotational  distributions 
were  first  smoothed  using  a  15-point  moving  average,  and 
then  extrapolated  to  7  =  0.5  and  J  —  480.5.  The  coefficients 


M(i  7  FIk*  rotational  distribution  ( unsmoothed  )  of  the  Hal  ptoduct  m 
v  K  produced  by  the  Ha  •  HI  beam  gav  reaction 


J  Chem  Phys  .  Vol  85  No  2  15  July  1986 


Noda  etal. :  Impact  parameter  dependence  of  Ba-t-  HI  reaction 


861 


FIG.  8  (a)  The  specific  opacity  function  for  ihe  Ba  •  HI 
•  Bal(v  =  8)  +  H  system,  assuming  there  is  no  velocity  dependence  for 
the  reaction,  (b)  The  experimental  (dot)  and  calculated  (solid  line)  rota¬ 
tional  distributions  for  the  reaction  studied  using  the  specific  opacity  func¬ 
tion  shown  in  (a). 


FIG.  9.(a)  The  specific  opacity  function  for  the  Ba  *  HI 
-•  Bal(r  -  8)  J-  H  system,  assuming  that  the  reaction  probability  varies 
asr,,./  (b)  The  experimental  (dot)  and  calculated  (solid  line)  rotational 

distributions  for  the  reaction  studied  using  the  specific  opacity  function 
shown  in  (a) 


{a, }  were  computed  by  the  least-squares  method  with  the 
constraint  that  the  coefficients  must  be  positives'  Since  this 
procedure  is  extremely  sensitive  to  the  “noise,”  the  resulting 
coefficients  are  smoothed  using  five-point  moving  average. 

The  "best-fit"  specific  opacity  function  obtained  in  this 
manner  has  a  peak  near  2.6  A  with  a  width  of  1.0  A 
(FWHM)  as  shown  in  Fig.  8(a).  The  smoothed  observed 
and  calculated  rotational  distributions  are  also  compared  in 
Fig.  8(b).  The  specific  opacity  function  in  Fig.  8(a)  gives  a 
remarkably  good  fit.  Thus,  it  is  concluded  that  for  this  model 
only  a  limited  range  of  impact  parameters  contributes  to  the 
formation  of  Bal(v  =  8). 

Very  often  the  overall  opacity  function  is  stylized  as  a 
step  function  or  a  smoothly  decreasing  function  of  the  im¬ 
pact  parameter  (see,  e.g.,  Fig.  6.6  of  Ref.  2).  These  simple 
models  have  been  successful  in  describing  impact  parameter 
dependence  in  a  previous  experimental  study”  on  Rb 
-i-  HBr.  However,  the  experimentally  determined  opacity 
function  reported  here  is  specific  to  a  particular  product  vi¬ 
brational  state,  namely  v  =  8,  while  the  often-quoted  opa¬ 
city  fi  nctions  are  not  for  a  specific  product  vibrational  state, 
but  for  all  product  states,  and  do  not  necessarily  resemble  the 
specific  opacity  functions. 


C.  Specific  opacity  function:  Intrinsic  velocity 
dependence  of  reaction  probability 

The  simplest  model,  considered  above,  ignores  any  ve¬ 
locity  dependence  in  the  reaction  probability.  We  consider 
next  how  the  specific  opacity  function  is  modified  if  there  is 
an  intrinsic  velocity  dependence.  We  have  investigated  two 
cases,  one  in  which  the  reaction  probability  varies  as  t!roll/:!, 
the  other  as  u]'f  The  resulting  specific  opacity  functions  are 
presented  in  Fig.  9  for  the  former  case,  and  in  Fig.  10  for  the 
latter. 

The  shapes  of  these  specific  opacity  functions  are  not 
very  different  from  the  specific  opacity  function  with  no  ve¬ 
locity  dependence  discussed  in  the  previous  section.  The  sig¬ 
nificant  difference  is  that  the  specific  opacity  functions  shift 
their  peaks  toward  small  impact  parameter  as  the  velocity 
dependence  is  varied  from  t',,.,1  ’.  c".,  (no  dependence),  and 
iff.  This  can  be  explained  as  follows:  if  the  collisions  with 
the  slower  relative  velocity  are  heavily  weighted  by  rltl'  . 
then  the  collisions  with  the  larger  impact  parameters  must 
have  more  probability  to  yield  the  same  rotational  distribu¬ 
tion  as  produced  by  the  r'.,  weighting.  Whatever  the  exact 
form  of  the  mtimsic  velocity  dependence  is.  we  conclude 
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HG  lo.i  n )  The  specific  opacity  fuiK Ill'll  for  fhe  Ba  -+  HI 
.  Balfv  Si  •  H  system.  assuming  that  the  reaction  probability  varies 
as  i  ,,,  (b)  The  experimental  (dot)  and  calculated  (solid  line)  rotational 

distributions  lor  (he  reaction  studied  using  the  specific  opacity  function 
shown  in  (a) 


that  this  type  of  velocity  dependence  has  little  effect  on  the 
specific  opacity  function  for  the  reaction  studied  here. 

0.  Specific  opacity  function:  Transmission  velocity 
dependence  of  reaction  probability 

As  a  last  model,  we  consider  the  situation  where  the 
reaction  channel  has  an  activation  energy,  i.e.,  only  colli¬ 
sions  with  relative  velocities  greater  than  some  threshold 
value  lead  to  reaction  (see  Fig.  11).  The  procedure  to  calcu¬ 
late  the  probability  for  a  trajectory  to  surmount  the  barrier  is 
based  on  the  “line-of-centers”  model,1 :  and  is  described  in 
the  Appendix.  In  this  model,  we  only  consider  the  activation 
energy  barrier  for  the  system,  and  the  energetic  constraint  on 
the  relative  velocity  distribution  is  not  included. 

The  Ba  -+  HI  reaction  is  only  slightly  exothermic.  The 
excess  energy  available  to  the  Bal  reaction  product  is  given 
b> 

£<BaI>  =  Z>i;  (  Bal)  £?”  ( HI )  +  £(HI) 

-  E{  H )  ‘  (13) 

where  /.  denotes  an  internal  energy,  D  the  bond  energy, 
and  I.  „  the  translational  energy  of  the  reactants  (double 
prime )  or  the  products  t  single  prune )  for  the  reaction  un- 
det  study.  A.  (  Ba )  0.  A ■  III  i  k  l .  the  average  rotational 


E  tfant  (keal/mol) 

0  5  10  15  20 


FIG.  II.  The  velocity  distribution  used  for  the  analysis  with  the  energy 
barrier  set  equal  to  (a)  Okcal/mol,  (b)  0.5  keal/mol,  and  (e)  2.0kcal/mol. 
Since  the  height  of  the  centrifugal  oarrier  varies  with  the  impact  parameter, 
the  cutoff  velocity  changes  as  a  function  of  b  A  value  of  b  --  3  A  is  used  in 
preparing  this  figure 


FIG  12. (a)  The  specific  opacity  function  for  the  Ba  -  HI 
.  Haiti  X)  •  II  si  stem,  assuming  that  Ihereexisis.in. iitoalion  energy 
of  US  Teal,  mol.  ibi  The  experimental  (doll  amt  calculated  i  solid  line) 
rotational  distributions  fi  :  the  reaction  studied  using  I  he  specific  opacity 
function  show n  m  i  a  i 
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This  tendency  has  already  been  noted  in  the  previous  sec¬ 
tion,  but  the  effect  of  the  transmission  velocity  dependence  is 
much  more  pronounced  than  the  intrinsic  velocity  depen¬ 
dence,  mainly  because  all  collisions  with  translational  ener¬ 
gy  less  than  the  activation  energy  barrier  are  unable  to  pro¬ 
duce  Bal(  v  =  8).  Nevertheless,  the  specific  opacity  function 
does  show  a  central  peak  near  2.6  A  suggesting  that  a  small 
range  of  impact  parameters  are  particularly  effective  in 
forming  the  Bal  product  in  v  =  8. 

Mims,  Lin.  and  Herrrr*  have  studied  this  reaction  by 
measuring  the  Bal  angular  distribution  produced  under 
crossed-beam  conditions.  Their  analysis  indicated  there  is  a 
translational  energy  barrier  of  2.5+  1.0  kcal/mol.  Since 
their  analysis  did  not  take  into  account  the  product  vibra¬ 
tional  and  rotational  excitation,  we  reanalyzed  their  angular 
distribution  data.  The  results  obtained  indicated  that  the  an¬ 
gular  distribution  is  fairly  insensitive  to  the  translational  en¬ 
ergy  barrier  height  and  a  moderate  agreement  can  also  be 
obtained  by  taking  the  Ba  +  HI  barrier  height  to  be  zero. 
Siegel  and  Schultz’''  studied  the  Ba  +  HC1  and  Ba  +  HBr 
reactions  at  different  collision  energies  and  suggested  that 
there  is  little  or  no  activation  energy  in  these  analogous  sys¬ 
tems. 

Because  of  the  limited  amount  of  information  on  the 
velocity  dependence,  we  cannot  determine  presently  the  spe¬ 
cific  opacity  function  for  this  reaction  system.  It  can  be  spe¬ 
culated  that  the  intrinsic  velocity  dependence  has  little  effect 
on  the  specific  opacity  function,  while  the  possible  existence 
of  an  activation  energy  barrier  and  the  uncertainty  in  the 
exoergicity  of  the  reaction  could  greatly  alter  the  form  of  the 
derived  specific  opacity  function.  Clearly  further  experimen¬ 
tal  work  is  needed  in  order  to  establish  the  reaction  cross 
section  as  a  function  of  the  collision  energy. 


energy./);  (HI)  70.4  kcal/mol.  M  and  D\(  Bal )  has  been 
variously  reported  as  71  4  -  1 .0  kcal/mol  from  a  mass  spec- 
trometric  thermochemical  study.  1  72.4  kcal/mol  from  a 
chemiluminescence  study.  '  and  less  than  77.4  kcal/mol 
from  a  study  oi  excited  state  predissociation.  ''  We  take  the 
\  alue  D  i  Bal )  72.4  kcal/mol  but  with  the  understanding 

that  the  dissociation  energy  is  still  not  established.  With  this 
value 

A.'(Bal)  3.1  kcal/mol  ■>  A'lr.ln„.  (14) 

However,  about  3.4  kcal/mol  is  required  to  populate  the 
r  S  vibrational  level.  While  we  anticipate  that  is 

much  less  than  because  translational  energy  of  the 

reactants  appears  primarily  as  rotational  energy  of  the  pro¬ 
ducts.  it  appears  nevertheless  that  the  Ball  v  =  8)  channel 
might  have  a  small  energy  deficit  w  hich  might  be  made  up  by 
the  translational  energy  of  the  reactants.  Again  this  suggests 
the  possibility  of  a  small  energy  barrier  in  the  reaction  chan¬ 
nel  This  is  expected  to  affect  the  opacity  function  in  a  way 
similar  to  how  an  activation  energy  barrier  does. 

Figures  12  and  13  show  the  specific  opacity  functions 
for  the  cases  where  the  energy  barrier  is  equal  to  0.5  and  2.0 
kcal  mol.  respectively  .  As  the  energy  barrier  becomes  high¬ 
er.  more  of  the  slower  velocities  are  rejected  causing  the  reac¬ 
tion  probability  to  increase  at  smaller  impact  parameters. 
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APPENDIX 

In  order  for  a  trajectory  to  reach  a  transition  stale,  suffi¬ 
cient  translational  energy  must  be  provided  to  surmount  the 
energy  barrier.  The  latter  consists  of  an  activation  energy  . 
E,„  and  a  centrifugal  energy  ,  E<  ,  i.e., 

/•; . .  /•/„  .  /  ;, .  (An 

where  equals  //  c;j1/2.  Hie  centrifugal  energy  is  ex¬ 

pressed  using  Eq  (  8  )  as 
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where  we  have  assumed  that  the  maximum  of  the  energy 
barrier  lies  at  d 

From  Eqs.  (Al)  and  (A2),  we  obtain 

>£,,(1  -  b2/d2)  (A3) 

Thus,  because  of  the  hne-of-eenters  functionality, 1  ’  the  min¬ 
imum  translational  energy  needed  to  reach  a  transition  state 
becomes  a  function  of  the  impact  parameter. 

In  the  actual  computation,  the  value  of  d  is  set  to  6  A, 
and  the  velocity  dependence  around  the  cutoff  velocity  is 
replaced  by  a  gradually  decreasing  function  rather  than  one 
hav  mg  a  sharp  cutoff. 
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